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Abstract. 

Observations of interstellar extinction and polarization indicate that the interstellar 
medium consists of aligned non-spherical dust grains which show variation in the inter- 
stellar extinction curve for wavelengths ranging from NIR to UV. To model the extinction 
and polarization, one cannot use the conventional Mie theory which assumes the grains as 
solid spheres. We have used a T-matrix based method for computing the extinction efficien- 
cies of spheroidal silicate and graphite grains of different shapes (axial ratios) and sizes and 
used these efficiencies to evaluate the interstellar extinction curve in the wavelength range 
3.4 - O.lyum. A best fit linear combination of silicate and graphite grains of not very large 
axial ratio, fits the observed extinction curve reasonably well. 

We calculate the volume extinction factor V^, which is an important parameter from the 
point of view of the cosmic abundance, for the spheroidal grain models that reproduce the 
interstellar extinction curve. We find that the shape of the grains do not affect the volume 
extinction factor 

Finally we have also studied the extinction and linear polarization efficiencies for 
aligned spheroids. The results show that the shape of grains affects the linear polarization 
efficiencies considerably for various orientation angles of the spheroids. 
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1. Introduction 

The exact solution to Maxwell's equations to calculate absorption, scattering and extinction of 
electromagnetic waves by homogeneous isotropic spheres is provided by Mie theory (Van de 
Hulst, 1981). Spherical grain models have been used by many authors to evaluate the interstellar 
extinction curve (e.g. Mathis et. al. 1977, Draine and Lee 1984). However, the spherical grains 
Send offprint requests to: Ranjan Gupta, e-mail: ragSiucaa . ernet . in 
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can not explain the observed interstellar polarization tiiat accompanies interstellar extinction. Mie 
theory would also not be strictly applicable if the particles have anisotropic optical properties, as 
for example graphite. Using the method of separation of variables Asano and Yamamoto (1975) 
have given the analytical solutions for homogeneous isotropic spheroidal particles. The method 
of separation of variables is based on expanding the incident, internal and scattered fields in vec- 
tor spheroidal functions. For spheroids with large sizes, large refractive indices or absorbing par- 
ticles, the system of linear equations becomes large and it becomes a difficult mathematical and 
numerical problem. Recently, Voshchinnikov and Farafonov (1993) have improved the separa- 
tion of variables method which is applicable to large aspect ratios. Several numerical techniques 
are now being in practice to obtain scattering and absorption cross sections for nonspherical 
particles (e.g. Bohren and Hufl'man 1983; Mishchenko et. al., 2000). The T-Matrix technique is 
based on expanding the incident field in vector spherical wave functions (VSWFs) regular at the 
origin and expanding the scattered field outside a circumscribing sphere of a scatterer in VSWFs 
regular at far field region. The T-Matrix transforms the expansion coefficients of the incident 
field into those of the scattered field and can be used to compute any scattering characteristic of 
a nonspherical particle (for details on T-Matrix method see Mishchenko et. al., 2002). 

One can also use the discrete dipole approximation (DDA see viz. Draine 1988) to study the 
extinction properties of the nonspherical inhomogeneous particles (see Borghese et. al. 2003). We 
have used the DDA to study the extinction properties for porous and composite particles ( Vaidya 
et. al. 1997, 1999, 2001; Paper 1, 11 and 111 hereafter). However, the DDA method required con- 
siderable computer time and memory. The T-Matrix code on the other hand (Mishchenko et. al., 

2002) runs much faster and the results obtained can be tuned with much ease since the input 
parameters to the code can be adjusted and re-run in a short time. However, the T-Matrix code 
in its present form cannot be used for studying inhomogeneous (e.g. porous, fluffy, composite) 
particles (Mishchenko et. al., 2002). 

In the section 2 we present the spheroidal grain models and describe the T-Matrix method; 
in section 3 we give the results of our calculations; i.e. extinction curves as functions of grain 
size, shape and wavelengths, the model interstellar extinction curves, and comparison of these 
model curves with the observed interstellar extinction curve (Savage and Mathis 1979, Whittet 

2003) and the influence of particle shape on the determination of interstellar extinction and in the 
section 4 we provide the conclusions of the present study. 

2. T matrix method and Spheroidal grain models 

The T-Matrix method was first introduced by Waterman (1965) for studying electromagnetic 
scattering by single, homogeneous nonspherical particles. The standard scheme for comput- 
ing the T-Matrix for single homogeneous scatterers in the particle reference frame is called 
the extended boundary condition method (EBCM) and is based on the vector Huygens prin- 
ciple (Waterman 1971). The general problem is to find the field scattered by an object bounded 
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by a closed surface. The Huygens principle establishes the relationship between the incident 
field, the total external field (i.e. the sum of the incident and the scattered fields), and the sur- 
face field. Technically the incident and the scattered waves are expanded in regular and outgo- 
ing vector spherical wave functions (VSWFs); the concept of expanding the incident and the 
scattered waves in VSWFs and relating these expansions by means of T-Matrix has proved 
very useful. The great computational efficiency of the T-Matrix approach to study electro- 
magnetic scattering by nonspherical particles with various shapes and sizes has found numer- 
ous applications. Several FORTRtVN T-Matrix codes are available on the world wide web at 
://www.giss.nasa.gov/~cnium| The codes compute the complete set of scattering character- 
istics, the optical cross sections, expansion coefficients, and the scattering matrix for randomly 
oriented and aligned particles. 

In the present study the axial ratios (ratio of horizontal to rotational axes hereafter called as 
AR) of the spheroidal grains ranging from 0.5 to 2.0 have been taken. Please note that for prolates 
we have AR < 1 and for oblates we have AR > 1 . The aspect ratio e is the ratio of the larger and 
smaller axes as defined by Mishchenko et. al.. (2002) and thus both a prolate with AR=0.5 and 
an oblate with AR=2.0 will have the aspect ratio e-2. For spheres AR=1.0 and e=l. 

The observed interstellar polarization requires that interstellar grains be aligned and non- 
spherical in shape. Hence, in addition to the calculations on random orientations (with orientation 
averaging) of the spheroidal grains we have also performed calculations for the aligned spheroids 
at several orientation angles. 

3. Results: 

3.1 Extinction properties of Spheroidal Grains 

One of the main objectives of the paper is to study the extinction efficiencies of the spheroidal 
silicate and graphite grains, obtained using the T-Matrix calculations and see the effects of shape 
of the grains i.e. spheres, oblates and prolates and various sizes which are typical to interstellar 
grain sizes i.e. 0.001 to about 0.500/i. These materials (i.e. silicate and graphite) have been the 
ingredients for most of the spherical grain models (see for example Mathis et. al. 1977, Draine 
and Lee 1984, Weingartner and Draine 2001). We present the results for the spheroidal grain 
models; i.e. the extinction efficiencies Qext as a function of wavelength in the spectral range (3.4 
to 0.1 ymi). 

Sihcates: 

The Figure 1 shows the extinction efficiency Qext for silicate grains for four grain sizes i.e. 
0.01, 0.05, 0.1 and 0.2 The y-axis in these plots are in the form of an extinction efficiency 
ratio Qext(Spheroid)/Qext(Sphere) to bring out the efltect of oblates and prolates as compared to 
spheres. 

It is seen from this figure that the grains with the same aspect ratio e-2 but having different 
axial ratios, AR=0.5 (Prolates) and AR=2.0 (Oblates), display similar extinction properties in the 
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Fig. 1. Extinction efficiencies (Qext(Spheroid)/Qext(Sphere)) versus wavelength for randomly ori- 
ented spheroidal silicate grains for four different grain sizes and five different axial ratios. 

entire wavelength range, i.e 3.4 - 0.1 jum. Results for the grains with other aspect ratio e =1.33, 
having an axial ratio AR=0.75 (Prelates) and AR=1.33 (Oblates) also show similar extinction 
property. These results on spheroids are consistent with the results reported by Mishchenko et. 
al. (2002). It is also seen that the extinction efficiency ratio Qext(Spheroid)/Qext(Sphere) for the 
silicate grains with aspect ratio e=1.33 (i.e. AR=0.75 and 1.33) is close to AR=1.0 (sphere), i.e. 
there is no appreciable variation in the extinction efficiency from spheres; whereas spheroidal 
grains with aspect ratio e=2 (i.e. AR=0.5 and 2.0), show considerable variation in the extinction 
from sphere i.e. the ratio Qext(Spheroid)/Qext(Sphere) deviates from 1.0. 
Graphites: 

Figures 2 shows extinction curves for spheroidal graphite grains. For small graphite grains 
(a=0.01 ju) the '2175A ' extinction feature is displayed by all the spheroidal grain models, but 
the peak is shifted for various axial ratios. In order to emphasize this aspect we show in Figure 
3 the extinction efficiency Qext for the spheroidal graphite grains for a=0.01 and 0.05 ii in the 
wavelength range 0.25 - 0.20 iim (i.e. 4-5 fim^^). These curves clearly show the shift in the 
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Fig. 2. Extinction efficiencies (Qext(Spheroid)/Qext(Sphere)) versus wavelengtii for randomly ori- 
ented spheroidal graphite grains for four different grain sizes and five different axial ratios. 



extinction peak with the shape of the grains (i.e. axial ratio). For a=0.01 the extinction peak 
at for spherical grains (AR=1 and Mie) shifts to for the spheroidal grains 

with AR=2.0. These curves also show the variation in the width of the feature. These results on 
the spheroidal grains indicate that the shape of the grains plays an important role in studying 
extinction properties and needs to be studied in more details. Voshchinnikov (1990) has also 
found the variation in the '2175A ' feature with the shape of the grain. In our earlier study on 
the porous grains (Vaidya et. al. 1997 & 1999 i.e. Paper I and II) we had found the shift in the 
the '2175A ' peak as well as the variation in the width of the feature for the graphite grains with 
porosity. Draine and Malhotra (1993) have found the shift in the central wavelength of the feature 
in the coagulated graphite grains but did not find any appreciable change in the width. It is to be 
noted from Figure 2 that for spheroidal graphite grains with aspect ratio e=1.33 (i.e. AR=0.75 
and 1.33), there is no appreciable variation in the extinction from spheres. It is also seen that 
for spheroidal graphite grains with a> O.OSjU, the extinction does not deviate much from spheres 
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Fig. 3. Extinction efficiencies Qext versus l/A for randomly oriented graphite grains of sizes 
a=0.01 and 0.05 jj. in the UV peak region for various axial ratios. The axial ration AR=1.00 
and Mie curves overlap each other as expected. The shift in the peak extinction is clearly seen 
for various axial ratios. 



in the optical and the UV spectral range. Please note that for prolate graphite with AR=0.5 in 
Figure 2 for a=0.01 and 0.05//, the T-Matrix calculation does not converge beyond 8 fim'^ . 
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Extinction and Polarization efficiency for Aligned Spheroids 

As mentioned earlier the existence of interstellar polarization requires that the grains must 
be nonspherical and aligned (see e.g Greenberg 1968, Wolff et. al. 1993, Voshchinnikov & 
Farafonov 1993). We have calculated the extinction and linear polarization for the aligned 
spheroidal silicate and graphite grains at several orientation angles. Figures 4(a-d) show Qext = 
[Qext(E) + Qext(H)]/2 and linear polarization efficiency IQpoil = Qext(E) - Qext(H) for aligned 
(at a fixed orientation angle /? = 45°) oblate spheroids for a = 0.1 //in the wavelength region 
3.4 - 0.3/(m; where Qext(E) and Qext(H) are extinction efficiency factors for the directions of the 
incident field vector Q(E) and perpendicular Q(H) to the axis of the spheroid. These results on 
the aligned oblates show that there is no appreciable variation in the extinction but the polariza- 
tion shows considerable variation with the axial ratio. Figs 4(e-h) show the extinction and linear 
polarization efficiency for a fixed axial ratio AR= 1.33 but with varying orientation angles /3 from 
to 90° for oblates. Here again the extinction does not show any appreciable change but the 
polarization does show a lot of variation. The aligned oblate silicate spheroids with AR =1.33 
produce maximum linear polarization between 2 and 3 yum ' , whereas the aligned graphite grains 
produce the maximum polarization between 1 and 2 /im '. These results on linear polarization 
are consistent with the results obtained for spheroids by Wolff et. al. (1993). Figs 5(a-h) show ex- 
tinction and Unear polarization efficiencies for the aligned prolate silicate and graphite spheroidal 
grains. The results in 5(a) and (b) on aligned prolates do not show appreciable variation with the 
axial ratio. However, the extinction varies with the orientation angles in 5(e) and (f) for aligned 
prolates. Note that the polarization efficiency of oblate silicates in Figures 4 (c & g) and prolate 
silicates in Figures 5 (c & g) show considerable differences but rest of the curves for oblates and 
prolates follow similar trends. 

3.2 Interstellar Extinction Curve 

We use the extinction efficiencies of the siUcate and graphite grains and the power law grain 
size distribution (i.e. n(a) ~ a"^-^, MRN size distribution (see Mathis, et. al., 1977 and Paper 
I, II & III) to reproduce the average observed interstellar extinction curve (Savage and Mathis 
1979, Whittet, 2003). We have evaluated the interstellar extinction curves for grain size distri- 
bution with the smaller size limit starting at a=0.001/i and the larger size limit up to a=1.0ju. 
Normalized interstellar extinction curve is generated for both silicates and graphites and then 
combined linearly as described in our Papers I, II and III, to obtain the best fit with the observed 
extinction curve using ax^ minimizing scheme. 

Apart from the various sizes tried for the lower and upper Hmits of the size distribution; the 
power law index of -3.5 was also changed from -2.5 to -4.5 in steps of 0.5 and the value of -3.5 
seems to be the best of the lot for obtaining best fits to the observations. Finally calculations were 
also repeated with different axial ratios and the best fit value ofx^ = 0.012837 was obtained with 
AR=1.00 i.e. spheres as was predicted by the MRN model. However, the spherical grains are not 
realistic, and besides, the observed interstellar polarization requires the interstellar grains to be 
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Fig. 4. Extinction Qext and linear polarization efficiencies versus XjA for aligned oblate spheroids 
(Silicates and Graphites with grain size a=0.1//). The panels (a), (b), (c) and (d) are for a fixed 
orientation (j6=45°) and various axial ratios AR of the aligned spheroids. The panels (e), (f), (g) 
and (h) are for a fixed axial ratio AR=1.33 but with various orientation angles yS. 

non-spherical. Thus we keep the oblate spheroids with AR=1.33, for which the best ^X)^ value 
is 0.013210, as the model for interstellar extinction. 

The Figure 6 shows the best fit minimized model curve along with the observed inter- 
stellar extinction curve. The top panel show the full curve in the wavelength region 3.4 to 0.1 
jum i.e. far UV to NIR. The bottom panel highlights the UV bump region. The best model curve 
thus consists of a size distribution a=0. 005-0.250 \i in steps of 0.005 // and with spheroids with 
AR=1.33. The observed points (shown as filled square dots) have been evaluated at wavelength 
steps of 50A by interpolation since the original data (Savage and Mathis 1979, Whittet, 2003) 
has much fewer wavelength points though it covers this large range. Similarly the T-Matrix model 
calculations were specially performed at this high wavelength resolution to provide smoother in- 
terstellar curves. Such calculations are feasible with T-Matrix due to its intrinsic property for fast 
computations. 
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Fig. 5. Extinction Qext and linear polarization efficiencies versus 1 //I for aligned prolate spheroids 
(Silicates and Graphites with grain size a=0.1/i). The panels (a), (b), (c) and (d) are for a fixed 
orientation (fi=45°) and various axial ratios AR of the aligned spheroids. The panels (e), (f), (g) 
and (h) are for a fixed axial ratio AR=1.33 but with various orientation angles /3. 

3.3 Volume Extinction factor and Abundance constraints 

An important parameter from the point of view of cosmic abundance is the volume extinction 
factor Vc which is defined as the ratio of total volume of the grains to the total extinction at a 
given wavelength i.e. 2 W Z CextiX) (see for example, Greenberg 1968, Vaidya et. al. 1984). Vc 
directly determines the amount of material required to render the level of extinction at a specific 
wavelength. Table 1 shows the volume extinction factors Vc at A = 0.55jum for the spheroidal 
grain models that reproduce the observed interstellar extinction (Figure 6). 

Results in Table 1 indicate that the variation in Vc due to shape (oblates or prolates) for both 
silicate and graphite grains is less than 10 % , hence the atomic-abundance constraints based on 
values of Vc for either spheres or spheroids of each material would be identical. These results 
showing eff"ect of shapes on volume extinction factor are consistent with the results obtained by 
Greenberg (1978). 
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Fig. 6. The top panel shows the comparison of observed interstellar extinction curve with the 
best fit model combination curve for spheroidal silicate and graphite grains (AR=1.33) in the 
wavelength region of 3.4 to 0.1 /iw using a grain size distribution of 0.005-0.250 jj. in steps of 
0.005 /i. This plot has a wavelength resolution of 50A for the model values and also for the 
observed values though the observed values are available at much lower resolution and at much 
fewer wavelength points. The lower panel highlights the U V bump region of the interstellar curve 
fitting. 

Table 1. Volume Extinction factor Vc for spheroidal grains. 
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4. Summary and Conclusions 

The extinction efficiencies for prolates, oblates and spheres are obtained using the T-Matrix cal- 
culations for silicate and graphite grains (of different sizes) in the wavelength range of 3.4 - 0.1 
/ZOT. The study brings out the effect of the shape of the grains on the extinction and Unear po- 
larization. Our results show that the shape of the grain plays an important role in detemnining 
the extinction properties of the grains and requires a more detailed study with other shapes and 
axial ratios. The spheroids with same aspect ratios (e) tend to show similar deviations in Qext 
values when compared with the Qext values of spheres and this phenomena has been noted by 
Mishchenko et. al. (2002). 

The Qext computations were used to model the interstellar extinction and then compared with 
the observed curve. Best fits are obtained for grain size distribution a=0.005-0.250ju in steps of 
O.OOSjU (with the power law exponent value of -3.5) with spheroidal grains having axial ratio of 
AR=1.33. 

The bare graphite/silicate grain model which we have used in the present study is consis- 
tent with observations of interstellar extinction. The volume extinction factors determined for 
the spheroidal grain models do not show appreciable variation from spherical models, so the 
abundance constraints for spheroidal grains and spheres would be identical. It should be noted 
the bare graphite/silicate grain model requires higher abundances of carbon and silicon than are 
recentiy available values for the interstellar medium (Mathis 2000). The real interstellar grains 
might be more complex (composite, porous; see e.g. Mathis 1996, Vaidya et. al. 2001). 

It is to be noted here that the T-Matrix method to study the optical properties of the spheroidal 
siUcate and graphite grains is a new approach and there is enough scope for varying the grain 
shape, size and composition for obtaining better model fits to the observations. Recentiy com- 
posite grains have been used by Moreno et. al. (2003) to model Comet Hale-Bopp dust grains. 
Petrova et. al. (2000) have used modified version of T-Matrix to study grains and have applied 
these results to model cometary grains. 
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